Properties of the intermediate type of 
gamma-ray bursts 



I. Horvath*, F. Ryde \ L.G. Balazs **, Z. Bagoly* and A. Meszaros § 

* Department of Physics, Bolyai Military University, Budapest, Box-12, H-1456, Hungary 
^Stockholm Observatory, AlbaNova, SE-106 91 Stockholm, Sweden 
**Konkoly Observatory, Budapest, Box-67, H-1525, Hungary 
* Laboratory for Information Technology, Eotvos University, Budapest, Pdzmdny P. s. 1/A„ H-1117, 

Hungary 

^Astronomical Institute of the Charles University, V Holesovickdch 2, CZ 180 00 Prague 8, Czech 

Republic 

Abstract. Gamma-ray bursts can be divided into three groups ("short", "intermediate", "long") with 
respect to their durations. The third type of gamma-ray bursts - as known - has the intermediate 
duration. We show that the intermediate group is the softest one. An anticorrelation between the 
hardness and the duration is found for this subclass in contrast to the short and long groups. 
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TWO-DIMENSIONAL GAUSSIAN FITS 

Simultaneously Mukherjee et al. [Q]] and Horvath [Hi found a third group of gamma-ray 
bursts (GRBs). Somewhat later several authors [HTH, d, ® also suggested the existence 
of the third ("intermediate") group as well. The physical existence of the third group is, 
however, still not convincingly proven. For example, Hakkila et al. IH] believe that the 
third group is only a deviation caused by a complicated instrumental effect, which can 
reduce the durations of some faint long bursts. Later Hakkila et al. tn\ published another 
paper which had different conclusions. 

Using Principal Component Analysis (PCA), Bagoly et al. [8] have shown that there 
are only two major quantities necessary (called the Principal Components; PCs) to char- 
acterize most of the properties of the bursts in the BATSE Catalog. Consequently, the 
problem of the choice of the relevant parameters describing GRBs is basically a two- 
dimensional problem. For the statistical analysis the choice of two independent parame- 
ters is enough; they may be, but are not necessarily, the two principal components. This 
means that only two parameters, relevantly chosen, should be enough for the classifica- 
tion and determination of the groups. Here we have chosen duration Tgo and hardness 
H32 — F3/F2 (F3 and F2 are the fluences) for these parameters. 

We can assume that the observed probability distribution of GRBs in this plane is a 
superposition of the distributions characterizing the different types of bursts present in 
the sample. Introducing the notations x = log^o and y = log #32 and using the law of 
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FIGURE 1. The 1956 GRBs in the {logTijojlog/^} plane. The lc ellipses of the three Gaussian distri- 
butions are also shown, which were obtained in the ML procedure. The different symbols (crosses, filled 
circles and open circles) mark bursts belonging to the short, intermediate and long classes, respectively. 



full probabilities we can write 

k 

p(x,y) = Y<p( x iy\ l }pi- (1) 
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In this equation p(x,y\l) is the conditional probability density assuming that a burst 
belongs to the l-th class. p\ is the probability for this class in the observed sample, where 
k is the number of classes. In order to decompose the observed probability distribution 
p(x,y) into the superposition of different classes we need the functional form of p(x,y\l). 
The probability distribution of the logarithm of durations can be well fitted by Gaussian 
distributions, if we restrict ourselves to the short and long GRBs, respectively [|9D. We 
assume the same for the y coordinate as well. Therefore it holds 
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The observational data from the Current BATSE GRB Catalog will be used. There are 
2702 GRBs, for 1956 of which both the hardnesses and durations are measured. 
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TABLE 1. Parameters of the Gaussian fits (k = 3) in the {log Tgo; log H32} plane. 



classes 


Pi 


a x = log T 90 


a.v = log #32 


o x 


<7y 


r (corr.coef.) 


1 


.245 


-.301 


.763 


.525 


.251 


.163 


2 


.109 


.637 


.269 


.474 


.344 


-.513 


3 


.646 


1.565 


.427 


.416 


.210 


-.034 



In order to find the unknown constants in Eq.© we use the Maximum Likelihood 
(ML) procedure of parameter estimation ifioll . One can define the Likelihood Function 
in the usual way, after fixing the value of k, in the form L = Elogp(jc ; ,y ; ) , where 
p(xi,yi) has the form given by Eq.®. Similarly, as it was done by Balazs et al. Q10|], 
the EM (Expectation and Maximization) algorithm is used to obtain the a x ,a y ,a x ,G y ,r 
and pi parameters at which L reaches its maximum value. We made the calculations for 
different values of k in order to see the improvement of L as we increase the number of 
parameters to be estimated. 

The confidence interval of the parameters came from 2(L max — Lo) = %m equation. 
Moving from k=2 to 3 the number of parameters m increases by 6 (from 1 1 to 17), and 
L,nax grows from 1193 to 1237, which means a very low probability of being a chance 
OCT 10 ). 

Moving from k=3 to 4, the improvement in Lmax is only 6 (from 1237 to 1243), which 
can happen by chance with a probability of 6.2%. Hence, the inclusion of the fourth 
class is not justified. Table 1. shows the parameters of the best fit with k=3. 



RESULTS 

The mathematical deconvolution of the p(x,y) joint probability density of the observed 
quantities into Gaussian components does not necessarily mean that the physics behind 
the classes obtained mathematically is different. It could well be possible that the true 
functional form of the distributions is not exactly Gaussian and that the algorithm of 
deconvolution formally inserts a third one only in order to get a satisfactory fit. One 
needs detailed investigations based on the physical (e.g. spectral) properties of the 
individual bursts to prove its astrophysical validity. 

Norris et al. f| 1 111 and Balazs et al. tlOfl found compelling evidence that there is 
a significant difference between the short and long GRBs. This might indicate that 
different types of engines are at work. The relationship of long GRBs to the massive 
collapsing objects is now also observationally well established, and the relation between 



the comoving and observed time scales is well understood |12L I13H . The short bursts 
can be identified as originating from neutron star (or black hole) mergers. Therefore 
the mathematical classification of GRBs into the short and long classes is physically 
justified. 

An important question that must be answered in this context is whether the interme- 
diate group of GRBs, obtained in the previous paragraph from the mathematical classi- 
fication, really represents a third type of burst physically different from both the short 
and the long ones. 



The classification into the short, intermediate and long classes is based mainly on the 
duration of the burst. From Table 1 one may infer that these three classes differ also in 
the hardnesses. The difference in the hardnesses between the short and long group is well 
known fbH] . According to these data the intermediate GRBs are the softest among the 
three classes. This different small mean hardness and also the different average duration 
suggest that the intermediate group should also be a different phenomenon, that is, both 
in hardness and in duration the third group differs from the other two. On the other 
hand, no significant correlation exists between the hardness and the duration within the 
short and the long classes. Thus, these two quantities may be taken as two independent 
variables, and the short and long groups are different in both these independent variables. 

In contrast, there is a strong anticorrelation between the hardness and the duration 
within the intermediate class. This is a surprising, new result, and because the hardness 
and the duration are not independent in the third group, one may simply say that only 
one significant physical quantity is responsible for the hardness and the duration within 
the intermediate group. Consequently, the situation is quite different here, because one 
needs two independent variables to describe the remaining two other groups. This is a 
strong constraint in modeling the third group. Hence, the question of the true nature of 
the physics in the intermediate group remains open, and needs further analysis. 

In this paper we have shown that statistically a third group of GRBs exists. Also 
statistically no further groups are needed to describe the {log Tgo; \0gH32} distribution of 
bursts. Finally, 1 1% of GRBs in the Current BATSE Catalog belong to the intermediate 
class. The memberships of the cataloged bursts are available on the internet jl5\ \. 
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